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lepton-flavor violation. 
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1 Introduction 



One of the main virtues of experiments at an e + e~ linear collider (LC) is the clean environ- 
ment allowing studies of the production and decay of new particles with low background. 
This not only enables precision measurements of particle properties, but also searches for 
very rare processes and small effects. An important example of this kind is lepton-flavor 
violation (LFV) as suggested by the experimental evidence for neutrino oscillations P3 |2] 
and expected particularly in SUSY models. Phenomenological investigations have indicated 
how tests of LFV at a high-energy LC could nicely complement searches for lepton-flavor 
violating rare decays such as \x — > ej. Previous work jSj has mainly focussed on slepton-pair 
production assuming two-generation slepton mixing. LFV in chargino-pair production was 
analyzed in |3] together with a detailed background analysis. Recently, it was shown in 
a rather model-independent analysis 5j by scanning over all possible soft SUSY-breaking 
terms consistent with the existing low-energy bounds that LFV in the above processes may 
be quite sizable, leading to final states with a pair of charged leptons of unequal flavor, 
missing energy, and possibly additional leptons and jets. 

In supersymmetric theories with heavy right-handed Majorana neutrinos, the seesaw mech- 
anism jHj can give rise to light neutrino masses at or below the sub-eV scale. Moreover, the 
massive neutrinos affect the renormalization group running of the slepton masses, which 
leads to mixing of different slepton flavors. In the present article we investigate the impli- 
cations of recent neutrino measurements on this mixing, extending our previous analysis of 
the radiative decays l a —>■ Ipj in the SUSY seesaw model [7j to the lepton-flavor violating 
processes e ± e~ — ► Ip^aXbXa involving slepton-pair production and subsequent decay. We 
again use the mSUGRA benchmark scenarios proposed in [S] for LC studies, concentrating 
on those which predict charged sleptons that are light enough to be pair-produced at the 
center-of-mass energy y/s = 500 GeV. Furthermore, we examine the most important back- 
ground processes and indicate the most promising lepton-flavor violating channels. The 
fact that in the benchmark models considered LFV occurs only in the left-handed slepton 
sector makes the separation of signal and background more difficult than in models with 
LFV in the right-handed sector. Therefore, our choice of models provides good study cases 
for developing search strategies. Finally, we work out the correlations between LFV in the 
high-energy e ± e~ collisions and the radiative lepton decays, and show that these correlations 
are relatively weakly affected by the uncertainties in the neutrino data, but very sensitive 
to the mSUGRA parameters. Consequently, they could play an important role in probing 
the class of models of LFV studied here. 
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The paper is organized as follows. In section 2 we briefly outline the SUSY seesaw mech- 
anism and discuss the slepton and neutralino mass matrices. Section 3 summarizes our 
analytic results on the helicity amplitudes for slepton-pair production, e ± e~ — > ifl^, and 
slepton decay, if — > l^Xa-i an d provides the cross-section formula for the complete 2 — > 4 
processes e ± e~ — > IplaXbXa- The fundamental parameters of the mSUGRA benchmark sce- 
narios as well as the neutrino data used in our numerical analysis are specified in section 4. 
In this section we also present our predictions on the signal cross-sections, while the esti- 
mates of the background and the effects of energy and angular cuts are given in section 5. 
The prospects for LC searches anticipated on the basis of our analysis are summarized in 
section 6. 

2 Lepton-flavor violation in the charged slepton sector 
2.1 Supersymmetric seesaw mechanism 

In supersymmetric models the existence of right-handed neutrino singlet fields with Majo- 
rana mass terms and Yukawa interactions gives rise to the well-known seesaw mechanism, 
provided the Majorana mass scale Mr is much higher than the scale v = 174 GeV of 
electroweak symmetry breaking. At energies below Mr the relevant term in the effective 
superpotential is given by jUj 

W e Jf = 1 -{Y U L ■ H 2 fM-\Y v L ■ H 2 ), (1) 

where L denotes the left-handed lepton doublets, H 2 the Higgs doublet with hypercharge 
+2 j an d Y u the matrix of neutrino Yukawa couplings. After electroweak symmetry breaking, 
the Yukawa couplings generate the Dirac mass matrix mo — Y U (H 2 ), (H 2 ) = v sin (3 being 

(H°) 

the H 2 vacuum expectation value with tan/? = jjgsj- This in turn leads to the Majorana 
mass matrix 

M u = m T D M- x m D = Yj M~ l Y u {y sin/?) 2 (2) 
for the light neutrinos which is diagonalized by the unitary MNS matrix U : 

U T M V U = diag(m 1; m 2 , m 3 ). (3) 

The matrix U and the mass eigenvalues rrii are constrained by neutrino data as discussed, 
e.g., in H|. 
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2.2 Renormalization group evolution of the slepton masses 

The heavy neutrino mass eigenstates contribute to the renormalization group running of 
the slepton mass matrices thereby inducing flavor non-diagonal terms which are responsible 
for the lepton-flavor violating processes described later. At the unification scale Mx we 
assume the mSUGRA universality conditions 

m\ L = mgl, mf R = m 2 l, A = A Y h (4) 

where mo is the common scalar mass and A Q the common trilinear coupling. At lower scales 
of order of the SUSY threshold, the mass squared matrix of the charged sleptons has the 
form 

„ / mf m~ ' \ 
m\ = **• ^ , (5) 
\ mf mf I 

\ Ilr Ir ' 

where mf , mf and mf are 3x3 matrices in flavor space, mf and mf being hermitian. 

'l i-r Ilr I'L Ir 

The respective matrix elements are given by 

(m-Jij = {m 2 L )ij + 6ij (mf. + m| cos 2/5 ^-^ + sin 2 6*^/^ (6) 
( m f R )ij = (rn 2 R ) ij + Sij (ml -m 2 z cos 2(3 sin 2 6 W ) (7) 
{m 2 LR )ij = Aijv cos (3 - 5ijm k fi tan (3, (8) 

9w being the weak-mixing angle, and fi the SUSY Higgs-mixing parameter. The contribu- 
tions to the first terms on the r.h.s. of © - © ar e indicated below: 

m 2 L = mil + (5m 2 L ) MSSM + 8m 2 L (9) 
m\ = mil + (5m 2 R ) MSSM + 5m 2 R (10) 
A = AoYi + SA MSaM + SA. (11) 

Here, (Sm 2 L r ) M ssm and 5A MSSM denote the usual MSSM renormalization group corrections 
[THj which are flavor- diagonal. In addition, the right-handed neutrinos radiatively induce 
the flavor off-diagonal terms 5m ^ and 5 A. In leading logarithmic approximation 1 , one 
gets [IH] 

5m 2 L = -^(3ml + A 2 )(Y}LY u ) (12) 



lr The exact one-loop renormalization group equations may give somewhat different results as pointed out 
in 1 1 1 1 and studied in more detail for lepton-flavor violating processes of interest in a forthcoming article 
|12j . For consistency with the results on radiative decays, derived in |Jj and used in the present analysis, 
we restrict ourselves here to the approximation l|12|) - (|14|l . 
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5m 2 R = (13) 

OA 

SA = -^JLOrrfLYj (14) 



' M 



with 

/",.-, = In i 

Mi 

Mi, i — 1, 2, 3 being the eigenvalues of the Majorana mass matrix M which may be chosen 
diagonal. Using (J2J) and (J3J), the Yukawa matrix Y v can be parametrized as follows [0]: 

Y v = — ^-diag ( JmI, Jm~ 2 , Jlfa) R diag (VmT, v 7 ^ v 7 ^) ^ (16) 



f sin (3 

where R is an undetermined complex orthogonal matrix. For substitution in (fT2*j) and (pPfj) 
the Yukawa couplings have to be evolved from the low-energy scale taken to be Mz to the 
Majorana scale Mr and further to the GUT scale Mx- 

Following [7j and many previous studies, we assume degenerate Majorana masses Mj.,2,3 = 
Mr and take R to be real. In this case, R drops out and the product YjLY^ is simply given 
by 

YjLY u = YX In i^f] = -^r^U ■ diagK, m 2 , m 3 ) • 17+ In (^*) . (17) 
V .Mr / tr sin p V Mr J 

A more general investigation allowing for non-degenerate Majorana masses is in progress 

and will be reported elsewhere [Tl|. Here, we stick to ()17j) since we want to investigate the 

correlations between the high-energy processes, which are the main subject of the present 

paper, and the low-energy rare decays studied in [7j. There, one can find some preliminary 

qualitative remarks on the effects of a complex matrix R. 

From (jSJ) the mass eigenvalues m ; ~ , i = 1, 6 of the sleptons are obtained via diagonaliza- 
tion by a 6 x 6 unitary matrix Up. 

,2 \ / 1 o\ 



U-m-Uj = diag(m| , ... 



,m- k ,...,m~ hi 



For definiteness, the eigenvalues are ordered such that the masses increase from li to m 
other words, in the absense of LFV l\ = f±, I2 = fii, h — t%. The corresponding mass 
eigenstates are then expressed in terms of the gauge eigenstates by 

k= {Ul)a.iha+ (Uf)(a+3)ilBa, i = 1, . .. , 6] a = e, fi, T. (19) 

2.3 Neutralino mass matrix 

In slepton production and decay neutralinos play an important role. Therefore, we want 
to clarify our notation for later use. The physical neutralinos are mixtures of gauginos and 
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higgsinos. In the gauge-eigenstate basis ip° = (B, W°, H®, H^) 7 , the neutralino mass term 
reads [T3] 

(20) 



(^°) T M^ U + h.c, 



with 



M,k 



Aft 






M 2 



-m z s w cp m z s w sp \ 
m z c w Cj3 -m z c w sp 
-p 








-m z s w c p m z c w Cf} 
\ m z s w sp -m z c w s l 3 —\i 

Mi being the U(l) and M 2 the SU(2) gaugino mass. The abbreviations s q 

= cos introduced above will be used throughout the paper. 
The neutralino mass eigenstates are given by 



Xa = N ab ipl a,6=l,...,4, 



where iV is the unitary matrix diagonalizing 

N^M^N' 1 = diag (m^o, 



(21) 



sin</> and 



(22) 



(23) 



and niyo, a 



neutralinos are composed of the Weyl spinors (j22|) as follows 



4 are the mass eigenvalues. Finally, the Majorana spinors for the 



(24) 



X a 



V° 
A a 

V° 
An 



3 Amplitudes and cross-sections 

The flavor mixing in the slepton sector induced by the heavy neutrinos as outlined in 
section 2.2 gives rise to the lepton-flavor violating processes e ± e~ — > lfl~ — > iflaXbXai 
i,j = 1,...,6. More specifically, in these processes LFV is caused by the slepton mixing 
matrix Uj in (|19|). which enters both the slepton production and the decay vertices. As 
a consequence, factorization in production cross-section times branching ratios is not ap- 
propriate. One rather has to coherently sum over the intermediate slepton states. In the 
following, we summarize our analytical results on the amplitudes and cross-sections for the 
above processes. The detailed calculations have been carried out in JE] and [T7j . 
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3.1 The lepton- flavor violating vertices 

The lowest-order Feynman diagrams for the processes under consideration are shown in 
Fig. and Fig. |2] with the particle four-momenta being defined in brackets. While the 
jli'lT vertex appearing only in Fig. Q is flavor-diagonal, LFV occurs in both the Zljl~ and 
/~/~Xa vertices described by 

-ZijiPs-PiT (25) 



with 



and 



with 



Cw-S\V 

1 3 

z ij = -n J2( u i)^(U~) ja + s^Sij (26) 

a=l 



- V2ie (A ia>a P L + B ia>a P R ) (27) 



in N* N* tw ~\~ N* 

A ia ,a = ~ ( U i)i(a+3) ~ „ ~ (#f)ia (28) 

ZmwSwCp Zsw 

Bia,a — —^-(Uj)i( a+ 3) + ~ " — ^—(Urjia, (29) 

cw zmwswcp 
respectively. Here, Pl,r = |(1 T75) are the left-right projectors, while the matrix elements 
N a b are defined in (j22J> and U{ a in (JTHJ). The lplj~Xb vertex is obtained from (J27|) after 
replacing v4 iaja by and B ia>a by A*^. 

3.2 Helicity amplitudes for slepton production and decay 

Although we will not consider e ± -beam polarization in the subsequent numerical analysis, 
we nevertheless present here the helicity amplitudes for future studies of polarization effects. 
In comparison to the slepton and neutralino masses, and to the cms energy of a LC the 
masses of the initial and final state leptons can be neglected in the kinematics. In terms 
of the Mandelstam variables s = (pi + P2) 2 , t = (pi — P3) 2 , u = (p\ — P4) 2 the helicity 
amplitudes Mij(h e ~, h e +) (h = ±1) for the production process e + e~ — ► 1^1^ are then given 
by 



•1 



M ij {+,+) = -2ie 2 V^£m*a Ba ' cA * jl ' c 



M y (--) = -2*e 2 v^£ 



$ij , 1 Zij , x B iX)C B* x 



e 2 — 1 y 4 4 A* 

b W 2 Zi 3 ■ ll < c jl,c 



*%(-,+) = +2ieVfa-^p + ^ 2 3 % +E 7 !i M L£ l- (30) 

\ s S VK C W / S — m Z c=l 1 — m x2 

Analogously, the helicity amplitudes Mij(h e -, h e -) for e~e~ — > read: 



M ij(->~) = 2ie 2 y/sJ2m^o 1 ' c : ' ' c + (t -» it, j <-» j) 



i 



4 



Furthermore, for the decay — > Z~x° the helicity amplitudes M~(h r ) summed over the 
helicity of the neutralino are given by 



M"(+) = 2ieB* a>a ^—p- a 

Mr(-) = 2zeA*^ a ^-p~ a . (32) 

From these one can obtain the helicity amplitudes M+(/ij+) for Il~ — > by the substitu- 
tions B ioija -> A* A6 and A iQ!)a -> S^ >6 [IS]. 

3.3 Cross-sections 

Having the amplitudes at hand it is straightforward to derive the cross-sections for the 
complete 2 — > 4 processes e ± e~ — > Z^/ 7 — > ZflZaX&Xa- F° r the square of the amplitudes 
summed over all possible intermediate slepton states one finds 



|M| 2 = Y.{M ij Ml l ){MrM^){MfM^)C ik7 ^C j 



vr „ 7T 



with 



x (*(pg - m? ) + <y(p§ - m|)) (5(pl - mf.) + ^ - mf )) , (33) 



r 1 



- , . Am 2 , ' 

1 



mr ifc = -(m r .r r .+m rft r r J, Am lk = m h - m Ik . (34) 



8 



Here, we have used the narrow width approximation for the slepton propagators, which is 
well justified since the slepton widths iy are of order GeV or less, that is much smaller 
than the slepton masses of order 100 GeV. For the product of two slepton propagators, this 
approximation yields (see, e.g., 

tA ( -i ^ — 1 « C lk ^^ (dip 2 -mf) + S(p 2 -mfj). 

p z — m~ +imiTij \p z — m~ + / 2mTik v 1 k ' 

(35) 

Moreover, for such small slepton decay widths one can also neglect the interference between 
identical outgoing neutralinos, as we have checked numerically. 
Integrating over the slepton momenta squared p\ and p\ and using the definitions 

, 4x , . d 3 p 3 d 3 p A MijM^ 



da ijkl = (2tt) 5{p x +P2-P3- Va) 



(2tt) 3 2£ 3 (2tt) 3 2£ 4 2s 



dB, = (2n)%p A -p & -p & ) / ^f 6 - ^f 8 £Lm±M±* (36) 

' ' y J (2vr) 3 2£ 6 (2vr) 3 2£ 8 2mr j7 7 z 



one can express the differential cross-sections in the form 



da = - ]T d(TijkidB ik dBji. (37) 



3 h' h 

„2 — ™2 2 



4 i; 



As one expects on general grounds, for large mass differences, Arh 2 k mT ik , the factors 
Cik in (}33|) approach 8^- Consequently, the coherent sum in (|37jl reduces to an incoherent 
sum over products of production cross-sections times branching ratios: 

da = X)d r(e ± e- - lfl-)dBr(l; -> r a xl)dBr(lf -> J±x°). (38) 

Previous studies often used only generic slepton masses and widths. In the present analysis 
the slepton decay widths r z ~. are actually calculated for each of the mSUGRA benchmark 
scenarios, including the effects of LFV. This is important for consistent phenomenological 
studies. 



4 Numerical results 

4.1 SUSY and neutrino parameters 

Similarly as shown in |7| for the radiative decays l a — > ^7, LFV in e ± e" — > ^^XbXa 
depends very sensitively on the SUSY scenario and on the neutrino masses and mixings. 
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The choice of models and the neutrino data are already described in some detail in the 
article mentioned above. Therefore, it suffices here to briefly specify the input parameters 
used in the following numerical analysis. 

Again, we choose the mSUGRA benchmark scenarios proposed in jS] for linear collider 
studies. In these set of models the universal trilinear coupling parameter A Q is assumed 
to vanish and, hence, LVF only occurs in the left-handed slepton sector as can be seen 
from (|12 |) -(|14 |) . For this reason, we focus on those scenarios which lead to sufficiently light 
sleptons, so that the production of slepton pairs containing a large admixture of at least one 
left-handed li component is possible at y/s = 500 GeV. The corresponding values of the 
mSUGRA parameters are listed in Tab. ^ together with some characteristic predictions. 



Scenario 


rai/2/GeV 


m / GeV 


tan/? 


m-jGeV 


VGeV 


m^o / GeV 


B 


250 


100 


10 


208 


0.32 


98 


C 


400 


90 


10 


292 


0.22 


164 


G 


375 


120 


20 


292 


0.41 


154 


I 


350 


180 


35 


313 


1.03 


143 



Table 1: Parameters of the relevant mSUGRA benchmark scenarios (from 8J). The sign 
of \x is chosen to be positive and A is set to zero. Given are also the mass and total width 
of the heaviest charged slepton and the mass of the lightest neutralino. 

Most likely, at the time when a linear collider will be in operation, more precise mea- 
surements of the neutrino parameters will be available than today. In order to simulate 
the expected improvement, we take the central values of the mass squared differences 
Am 2 ,- = rrij — mf and mixing angles 9{j from a global fit to existing data JB] with er- 
rors that indicate the anticipated 90 % C.L. intervals of running and proposed experiments 
as further explained in 

tan 2 6 23 = lAOtlil tan 2 13 = 0.005^, tan 2 12 = 0.36l£?jj, (39) 
Am 2 2 = 3.30l8:| • 10 " 5 ey2 > Am 23 = 3.10lJ;g • 10~ 3 eV 2 . (40) 

The fit [TH] is not the latest one, but still valid within the present uncertainties. It is taken 

for consistency with the results derived in |7| and used in the following. 2 The CP-violating 

2 In section 4.2 we will remark on the slight changes in our results when using the central values of the 
neutrino parameters from one of the more recent fits instead of |18j . 
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Dirac phase 8 is very difficult to measure and hence we allow it to vary in the full range 



5 = - 2tt. (41) 

Furthermore, we consider both hierarchical and quasi-degenerate spectra of light neutri- 
nos. The case of neutrino masses that are too light for an absolute mass measurement, is 
described by 



m\ = - 0.03 eV, m 2 = \jm\ + Amf 2 , m 3 = \jm\ + Am| 3 . (42) 

For mi <C m 2 <C 777,3, the combination YjY u , entering via (J 17)) the renormalization group 
correction (|12j) to the slepton mass matrix, can be approximated by 

( YX ) a w y^ u * u h + \^^:;,) • (43) 

For the quasi-degenerate case we take 

777-1 = (0.3+o}g) eV, 777 2 f« 777i H — Am? 2 , 777 3 Pi 777i H —Ami. (44) 

2/;/| - 2mi 



Since in this case mi ;§> y Am{ 2 - ! y Am 23 , one can use the approximation 

* A h' + + ^^)) • (45 » 

The Majorana mass M# is kept as a free parameter and is only constrained by the require- 
ment that Y v should stay small enough for perturbation theory to hold. 



4.2 Predictions for e + e and e e collisions 

Using the numerical input specified in the preceding subsection, we have calculated the 
cross-sections for e + e~ — > — > fi + e~ + 2%° and r + + 2%°, as well as the cross-section 
for e~e~ — > Ijl^ — >■ /i~e~ + 2x°. In the SUSY scenarios under study R-parity is preserved, 
and xl is the lightest SUSY particle and therefore stable. Since Xi is n °t detected, the 
observable final states consist of a pair of charged leptons with different flavor plus missing 
energy (Jfi). Thus one is dealing with a rather simple signal. Note that in the models defined 
in Tab. ^the heavier neutralino states a = 2,3,4 do not contribute to these signals: 
either they are too heavy and hence not or only very rarely produced in the slepton decays, 
or they are too light in order to decay invisibly via Xa ~ y vv — > Xi ul/ - Other open channels 
such as fr lead to more complicated final states which are not considered here. As a rule of 
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thumb, for an integrated luminosity of 1000 fb -1 , the goal at an e + e _ LC such as TESLA, 
the cross-section should be larger than 10 -2 fb for a signal to become observable. 
Fig. El shows the cross-sections for [i + e~ + 2x1 an d T + p,~ + 2\\ at ^/s = 500 GeV as a 
function of the right-handed Majorana mass scale Mr in scenario B which predicts the 
lightest sleptons. We have assumed a very light neutrino spectrum and have scattered 
the neutrino parameters in the error intervals specified in ([39 )1 - ()42 )) according to Gaussian 
distributions, except ni\ and S for which we have taken flat distributions. Each choice 
corresponds to one particular point in Fig. |H1 

A few comments concerning the main features of this plot are in order. Similarly as the 
branching ratios Br(l a — > Ipj) [7], the cross-sections exhibit the typical proportionality to 
(YJY u ) 2 qc Mft. This behavior indicates that the signal cross-sections are dominated by 
Feynman graphs involving the lepton-flavor violating vertices in first order. The latter in 
turn implies that e + e~ — > r + / u~ + 2x ( j ) is dominantly a s-channel process and thus suppressed 
relative to e + e~ — > [i + e~ + 2\\ where the s- and t-channel contribute. 
Furthermore, in contrast to the branching ratios of the radiative decays, the cross-sections 
saturate for large Mr. This saturation is seen most clearly in the /i + e~ channel. It can 
be understood by realizing that for large Mr the mass differences of the sleptons with a 
dominant left-handed component become comparable to the corresponding slepton widths. 
As already mentioned, in this case the cross-sections can be approximated by the incoherent 
sum (|38p. Focussing on the e/x channel, the dominant contribution is determined by the 
cl-^l mixing angle or, more precisely, by 



" ~ K)n - K> ~ (Khy n - (rJnW 



where the second equality follows from (jUJ), ©, and (fT5|) if the lepton masses are neglected. 
Using one then sees that tan26* eM and thus also the cross-sections are independent of 
Mr. In the r + /i~ channel the saturation is less pronounced due to effects of the heavier r 
mass which dominate the denominator in the expression analogous to (|46|). except for very 
large values of Mr. The finite r mass is also responsible for a suppression of all channels 
with a t lepton in the final state relative to the \ie channel. 

Finally, the impact of the neutrino uncertainties is weaker in the r + /i~ than in the /i + e~ 
channel. This is not surprising, since the product of neutrino Yukawa couplings (YjY v j^ 
relevant for r + /i~ is mainly given by a single term depending on the large angle 623 and 
mass squared difference Am| 3 , whereas \ Y^Y,/j relevant for jj + e~ depends on two terms 



involving the small quantities #13 and Am 2 2 , respectively, as can be seen in (|4"3j) and (|45j) 
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for i ^ j. Given the uncertainties in # 13 , the two terms can become similar in magnitude 
and cancel depending on the unknown phase 5. The same holds for the re channel. Since 
this channel is unfavorable because of the r mass suppression pointed out above and the 
comparably large MSSM background, we will concentrate on the /xe and r/x final states. 
The corresponding cross-sections for the higher cms energy y/s = 800 GeV are plotted in 
Fig.HJ Here, the scenario I being unaccessible at 500 GeV is included. This figure illustrates 
the range of cross-sections which can be expected in the benchmark models specified in Tab. 
HJ Together with Fig. |3] it also indicates the energy-dependence of the signals. 
If the central values of the neutrino parameters from one of the latest fits are used 
instead of (39) and (40), the most likely predictions for cx(e + e~ — > fi + e~ +2x1) T ^ se by about 
30% as compared to the results shown on Fig. 3 and 4, while the spread of the predictions 
shrinks slightly. Furthermore, the cross sections for e + e _ — > t + /j,~ + 2x1 decrease by about 
20% with the spread remaining unchanged. 

We now turn to a particularly interesting subject, namely the correlations between the 
cross-sections a(e ± e~ — > lfl~ — > IplZXbXa) anc ^ ^he branching ratios Br(l a — > Ipl)- This 
correlation is illustrated in Fig.EJfor e + e~ — > /x + e~ + 2x1 & t = 500 GeV and /x — > ej. For 
Mr < 10 13 GeV, corresponding roughly to Br(fi — > ej) < 10 -13 [Zj, this correlation is so 
accurate that the neutrino uncertainties drop out almost completely. Only for larger values 
of Mr the correlation gets lost because of the saturation of the cross-sections explained above 
which sets in at different values of Mr depending on the precise values of the neutrino 
parameters. Fig. |H1 shows the analogous relation between a(e + e" — > t + + 2x\) and 
Br(r — > yrf) which is even stronger. Note that beyond the upper ends of the scatter plots 
in Fig. El the neutrino Yukawa couplings would become too large for perturbation theory to 
hold. 

These correlations can be used to estimate the LC cross-sections that are allowed by bounds 
on or measurements of the radiative decays. Taking scenario B as an example, one can read 
off from Fig. El that the present bound Br(fi —>■ ej) < 1.2 • 10~ n [20J implies a(e + e~ — > 
/x + e~ + 2x1) < 0.3 fb, while a measurement of Br(fi — ► ej) ~ 10 -14 by the new experiment 
at PSI [21] would predict a(e + e~ —>■ /x + e~ + 2x1) ~ 3 • 10~ 4 fb. In other words, if /x — > c-f 
will not be detected at PSI, one does not expect an observable /xe signal at a 500 GeV LC 
either. In model C, the above branching ratios imply cr(e + e~ — * /i + e~ + 2xi) < 1 fb and 
pa 6 • 10 -3 fb, respectively. 

Analogously, from Fig. El one concludes that the present bound Br{r — > fij) < 6 • 10 -7 
[221 provides the weak constraints cr(e + e _ — > r + /x _ + 2x1) < 10 fb at y/s = 800 GeV for 
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model B, C, and G, and < 2 fb for model I. More importantly, also the sensitivity goal 
Br(r — > fi'-f) pa 10 -9 of future searches [23.] will not rule out sizable LC cross-sections, 
namely cr(e + e~ — > r + /i~ + 2x1) ~ 6 (0.2) fb in scenario C (B and G). In fact, in the 
Tfi channel LC experiments may reach farther than the r — ► /ry experiments planned in 
the future. At any rate, here we have a nice example for the complementarity of low and 
high-energy searches. 

In principle, one has such correlations also in channels which differ in flavor. This is exempli- 
fied in Fig.[7|for a(e + e~ — > r + fi~ + 2x1) and -Br(/i — > ery). However, because of the different 
flavor-violating couplings involved this correlation suffers considerably from uncertainties 
in the neutrino sector. Nevertheless, the experimental bound Br(fi — ► ej) < 1.2 ■ 10 -11 
yields a stronger constraint on model I than the one obtained from Fig. |Hl excluding 
a(r + fj,~ + 2x?) > few ■ 10~ 2 fb at ^ = 800 GeV already today. 

Finally, it is interesting to look at the relation between the two cross-sections <r(e + e~ — > 
r + / u _ + 2x1) an d cr(e + e~ — > /i + e~ + 2x?) themselves. This relation is displayed in Fig. |H] 
for 800 GeV and scenario C. Unfortunately, it is washed out by the uncertainties in the 
neutrino parameters. Nevertheless, a measurement in the fie channel would give a useful 
lower bound in the t/j channel. 

The use of the latest central values of the neutrino parameters from [19J instead of (39) and 
(40) does not change the correlations shown in Fig. 5 and 6 significantly. On the other 
hand, the lower boundaries of the correlations displayed in Fig. 7 and 8 are slightly shifted 
to the right. 

For completeness, we have also studied the prospects for e~e~ collisions. Fig. El shows 
<j(e - e - — > ji~e~ + 2x1) an d 0"(e~e - — > r~fi~ + 2x\) for yfs = 800 GeV and model C. 
The strong suppression of the t~/i~ final state results from the fact that at least two 
flavor-violating couplings are required because of the absence of the s-channel (see Fig. 
I2J). The r~e~ final state is suppressed for similar reasons as the r + e~ channel in e + e~ 
collisions. Thus, fi~e~ remains as the only promising channel. The correlation between 
<r(e~e~ — > jjTe~ + 2%?) and Br(n — > ej) is plotted in Fig. EH for the same scenario C. 
We see that the cross-section permitted by the present bound on /i — ► ej lies in the range 
of 1-10 fb which would make the detection rather easy. In this figure, we also show the 
prediction for a degenerate neutrino spectrum. Generally, both the LC cross- sect ions as 
well as the radiative branching ratios tend to be smaller for degenerate than for hierarchical 
neutrinos. This reduction comes from the suppression of the lepton-flavor violating terms in 
(|43|) by the neutrino mass mi which sets the absolute mass scale. However, the correlation 
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between cross-sections and branching ratios is found to be more or less independent of the 
absolute neutrino mass scale. The predictions shown in Fig. 9 would change only slightly, 
if the central values of the neutrino parameters from ^H] are used instead of (39) and (40), 
similarly as described above for e + e~ scattering. Again, the correlation depicted in Fig. 10 
remains unaffected. 

5 Background processes 

After having discussed the lepton-flavor violating signals leading to the channels lpl a +]fi we 
turn now to the lepton-flavor conserving standard model (SM) and SUSY background. For 
the dominant SM reactions we also illustrate the efficiency of angular and energy cuts in 
reducing this background. The SUSY background processes are calculated in the respective 
mSUGRA benchmark scenarios without cuts. We are going to assume that at the time 
these searches for LFV will be performed, SUSY has already been discovered and the 
lighter sparticle masses and main decay channels are known. In this case, one will be able 
to design very specific cuts optimized to each channel of interest. Under this assumption 
our estimates indicate that it should be possible to reduce the SM and SUSY background 
sufficiently and detect some of the signals considered, provided the Majorana scale Mr is 
not much smaller than 10 13 GeV. However, more conclusive statements on the feasibility 
to discriminate signal from background require detailed simulations which are beyond the 
scope of this paper. 

5.1 Background in e + e 

In the SM, the dominant background in the channels e + e~ — > l+l~ — > l^l~ + 2x < ( is produced 
by the following lepton-flavor conserving processes: 

Bl) e+ e - _> W+W- -> lp & l~v a (a ± /?), 

B2) e+e- -> W + e~v e -> l^upe~u e (p ^ 1), 

B3) e+e" -> t + t~ -> T + v T l-p a (a = 1, 2). 

These have been estimated with the help of the COMPHEP program package j21] , including 
the following kinematical cuts: 

i) beam-pipe cut: | cos(e ± , lp )\ < 0.985 (10°), 
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ii) lepton-energy cut: E min < E x < E max , 

iii) missing-energy cut: 2m^o <$ < \/s — 2E min . 

The angular cut i) partially eliminates the large contribution from t-channel photon ex- 
change (B2), and also small-angle W (Bl) and r (B3) production. Since the angular distri- 
butions of the decay leptons from the heavy sleptons are relatively flat, the signal rates are 
thereby reduced by less than 10%. Then, knowing the slepton spectrum and the LSP mass, 
one can impose the cuts (ii) and (iii) requiring the lepton energies and the total missing 
energy to lie in the intervals corresponding to the decay I — > lx%- Since one is dealing with 
a coherent process, the intermediate slepton flavor is not known. Therefore, E min (E max ) 
is defined as the minimum (maximum) of the decay energies kinematically allowed for any 
slepton flavor: 

/ m~ - m| \ ( m i ~ m l° \ 

Emm = min 2(£Ffl) ' Emax = max WF^i ' (47) 
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i+ii+iii, G 


11.6 




i+ii+iii, G 
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i+ii+iii, I 


10.8 




i+ii+iii, I 


24.2 


Final state 


v/i/GeV 


Cut 
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cr/fb 
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58.8 


800 


i 


22.5 






i+ii+iii, B 


11.1 




i+ii+iii, B 


6.0 






i+ii+iii, C 


3.1 




i+ii+iii, C 


3.5 






i+ii+iii, G 


3.5 




i+ii+iii, G 


3.9 






i+ii+iii, I 


3.6 




i+ii+iii, I 


4.8 



Table 2: Cross-sections for the SM background in e + e collisions including cuts as explained 
in the text. 

The SM background cross-sections that remain after these cuts are summarized in Tab. [21 
If one requires a signal to background ratio S/ \/B~ = 2 and assumes a typical signal cross- 
section of 0.1 fb, one can afford a background of about 2 fb. The additional suppression may 
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be achieved by applying the standard selectron selection cuts on the acoplanarity, lepton 
polar angle and missing transverse momentum. It has been shown that in this way the SM 
background to slepton-pair production can be reduced to about 2-3 fb at = 500 GeV, 
while the signal cross-section shrinks only by a factor 3. It should be noted that suppression 
of the W-pair background by right-handed electron polarization is not an option here, since 
in the mSUGRA scenarios considered LFV occurs only in the left-handed slepton sector, 
and thus right-handed polarization would also kill the signal. 
In the MSSM, the main background processes are 

B4) e+e~ — > v a v a , 

B5) e+e- -> l+!~, 

B6) e+e- -»■ xtXa , 

B7) e+e' -> x^e~z> e , 

with the sparticles decaying like v a — > v a Xi an d laXi ■> Xi l + v a and l+v a , l + - > ^Xi an d 
VaXii an d many other channels, depending on the given sparticle spectrum. 
The numerical estimates for the MSSM background listed in Tab. El are again obtained 
with the help of COMPHEP. Here, no cuts are included and the cross-sections for the 
individual production channels are added incoherently. We find that the MSSM background 
to \x+e~ + fi is very small, below 0.2 fb in all scenarios of Tab. Q except for model I and 
C, where it amounts to 0.4 and 5 fb, respectively, at y/s = 800 GeV. With 2-7 fb the 
r+fi~ background is considerably bigger. Also shown in Tab. El is the MSSM background to 
r+e~ + ft, which can contribute to the fi+e~ channel via the decay t+ — > \i+v il v T . If f\ and 
xf are very light, like in scenarios B and I, this background can be as large as 100 fb. The 
charginos originate mainly from selectron or v e pair production and decay almost exclusively 
into staus, Xi ~ * t~\V t . However, such events typically contain two neutrinos in addition 
to the two LSPs which are also present in the signal events. Thus, after r decay one has 
altogether six invisible particles instead of two, which should allow to discriminate the signal 
in n+e~ + If) also from this potentially dangerous MSSM background by cutting on various 
distributions. 

5.2 Background in e e 

In seesaw models the process e - e~ — > W~W~ — > lZl~Vpv a (a, (3^1) via neutrino exchange 
is highly suppressed due to the very small admixture of the heavy Majorana neutrinos in 
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Table 3: Cross-sections in fb for the MSSM background in e + e~ collisions at ^Js = 500 GeV 
and 800 GeV. The second column specifies the important production channels. In scenarios 
C and I the cross-sections are below 0.01 fb at y^s = 500 GeV and are therefore omitted. 
Otherwise, cross-sections below 0.01 fb are denoted by 0. Kinematically forbidden channels 
are marked by a hyphen. 
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i+ii+iii, I 


19.4 
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Table 4: Cross-sections for SM background processes in e e collisions including cuts as 
explained in the text. 

the light neutrino eigenstates [211]. The analogous argument holds for the SUSY version 
of this process, i.e. chargino-pair production via t-channel sneutrino exchange. The main 
background processes leading to the final state e~e~ — ► [i~e~ + $ are single W~ production 
(cf. B2), charged slepton-pair production (cf. B5), and single chargino production (cf. B7). 
The SM background e~ pTv^v^ originating mainly from single W~ production is estimated 
in Tab. E] including the cuts (i)-(iii) explained in the preceding subsection. As can be seen, 
the cross-sections are bigger and the cuts are less efficient than in the e + e~ case. Moreover, 
contrary to what we have found for e + e~ collisions, the situation at y/s = 800 GeV is 
somewhat worse than at 500 GeV due to the increase of single W~ production with energy. 
The total cross-sections for the MSSM background processes are summarized in Tab. El In 
all four scenarios the MSSM background to fi~e~ is below 1 fb, while the background 
to r~er in scenarios B and I is of the order of 100 fb. This is similar to our findings 
for e + e~ collisions. Despite the larger background, the signal to background ratio tends to 
be slightly more favorable in e~e~ collisions than in e + e~. For example, for scenario C and 
x/s = 800 GeV one finds roughly S/\^B(e~e~) ~ 2 S/y/B(e + e~), assuming an integrated 
luminosity of 250 and 1000 fb -1 , respectively. 

6 Conclusions 

We have studied the lepton flavor- violating processes e ± e~ — > J2ij iflT ~^ IplaXbXa a ^ a 
future linear e + e~ collider. As a theoretical framework we have chosen mSUGRA bench- 
mark scenarios furnished with the seesaw mechanism assuming a Majorana mass scale in 
the range 10 11 GeV < Mr < 10 15 GeV. Our analysis shows that for center-of-mass energies 
of 500-800 GeV, the signal cross-sections may be as large as 1-10 fb. The detailed predic- 



19 



Final state 


Production 


B 


C 


G 


I 


e~n~ +$ 




500 


800 


800 


500 


800 


800 




erer 


0.01 


0.01 





0.01 


0.2 


0.4 









0.01 


0.4 


- 


0.1 







Sum 


0.01 


0.02 


0.4 


0.01 


0.3 


0.4 


e~T- +$ 




500 


800 


800 


500 


800 


800 






140.8 


102.3 





0.2 


8.1 


60.0 









0.01 


0.4 




0.07 







Sum 


140.8 


102.3 


0.4 


0.2 


8.2 


60.0 



Table 5: Cross-sections in fb for MSSM background processes in e e collisions at y/s = 
500 GeV and 800 GeV as in Tab. El 

tions depend strongly on the SUSY parameters and on the neutrino masses and mixing. 
Consequently, they are also strongly affected by uncertainties in the neutrino data. 
On the other hand, the correlations between the high-energy cross-sections and the branch- 
ing ratios for the corresponding rare radiative decays are expected to be less influenced 
by experimental errors in neutrino parameters. We have therefore investigated these cor- 
relations very thoroughly and find that they are the stronger, the smaller the branching 
ratios for the radiative decays, that is the lower Mr. At Mr smaller than 10 13 GeV the 
neutrino uncertainties play no role at all. The present bounds on /i — > e'y and r — ■> /ry 
still allow sizable signals at a LC. If // — > e^y will not be observed at the new PSI exper- 
iment, the cross-section for e + e _ — > Y^ijlf^T ~^ fi + e~ + 2x\ is predicted not to exceed 
0.02 fb at yfs = 500 GeV. On the other hand, Brir — ► /ry) < 10~ 9 alone would still be 
compatible with a cross-section of order 1 fb for e + e~ — > IjlT ~ > t + + 2\i, whereas 
Br([i — > ej) < 10~ 14 would be in conflict with such a large cross-section. 
Finally, we have estimated the SM and MSSM background. With only a beam pipe cut 
of 10 degrees included, the SM background at 500 to 800 GeV amounts to 80 to 100 fb 
in the fi + e~ channel and to 20 to 60 fb in the r + /i~ channel. The corresponding MSSM 
background varies between 0.01 and 10 fb depending on the SUSY scenario, the final state 
and the beam energy. An exceptionally large MSSM background of 50 to 100 fb is found in 
scenarios B and I in the channel r + e~ which also contributes to the fi + e~~ channel through 
leptonic r decay. Since one can assume that the relevant sparticle spectra are known at 
the time such experiments can be done, it appears possible to sufficiently suppress the 
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MSSM background by carefully designed cuts on decay energies, missing energy, angular 
distributions and other quantities. However, conclusive feasibility studies require a detailed 
Monte Carlo simulation. 

The overall discovery potential for LFV at a LC is slightly increased by performing searches 
in e~e~ collisions in addition to e + e~ collisions. 
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Figure 1: Feynman diagrams for e + e — > — ► ^^XbXa- The arrows on scalar lines refer 
to the lepton number flow. 




Figure 2: Feynman diagrams for e e — > l i lpl a XbXa- The arrows on scalar lines refer 
to the lepton number flow. 
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gure 3: Cross-sections at y/s = 500 GeV for e + e — > /i + e + 2\\ (circles) and e + e" 
7i _ + 2^5 (triangles) in scenario B for the case of very light neutrinos. 
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Figure 5: Correlation of a(e + e — > fi + e + 2x1) a ^ = 500 GeV with Br{ji — > ej) 
scenario C (triangles) and B (circles) for the case of very light neutrinos. 
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scenario (from left to right) C, G (open circles), B and I for the case of very light neutrinos. 
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Figure 7: Correlation of <r(e + e — ► r + /i + 2x5) a ^ = 800 GeV with Br(fi — > ej) in 
scenario C (triangles) and I (circles) for the case of very light neutrinos. 
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Figure 8: Correlation of a(e + e — > r + ji + 2xi) with a(e + e — > ji + e + 2x°) at \fs = 800 
GeV in scenario C for the case of very light neutrinos. 
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gure 9: Cross-sections at yfs = 800 GeV for e e — ► \i e + 2x1 (circles) and e e" 
fi~ + 2xi (triangles) in scenario C for the case of very light neutrinos. 
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Figure 10: Correlation of <r(e~e~ — > yU~e~ + 2%°) at y/s = 800 GeV with Br(fi — > e7) 
in scenario C for the case of very light neutrinos (closed circles) and degenerate heavier 
neutrinos (open circles). 
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